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SUMEARY 


The object of this thesis is to determine the effect of small 
openings on the atress distribution in webs of beams subjected to jure 
bending and to devise a rapid means for predicting the increase in stress 


over that of the condition of no opening. 


Experimental evidence for a family of specimens which covered 
the usual structural practice hes been plotted in the RESULTS. A theor- 
etical solution, daveloped in the APPENDIX, is included end compared with 
the experimental avidence. From these plots the authors conclude thats 
the theoretical golution compares quite favorably with actual experiment; 
the stress concentration factor for the edge of the hole varies directly 
with the ratio of depth of opening to depth cf web and inversely with 
the ratio of length to depth of opening; and thut the stress concentration 
factor for the edge of the web varies directly with the ratio of depth 
of opening to depth of web, but the variation with the ratio of lenth 
to depth of opening is more complex showing a pronounced maximum in the 
vicinity of 1.5. The fact that the ratio of denth of opening to depth 
of web is of little importance until it reaches the value of about 1/3 


is considered excellent verification cf presently accepted practice, 


On the basis of the foregoing, it is considered highly desirable 
that further investigntion for specimens having a ratio of length to depth 
of opening in the neighborhood of 1.5 be mide with a view to confirmation. 


It is further considered desirable that a systematic survey cf a similar 


x 













family be made to determine the effectiveness of reinforcement at oj :n- 


ings and simultaneously cases of complex bending; 1.e., bending with shear. 
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INTRODUCTION 


In almost all structures cpenings of some sort are required, 
either for access of for passing piping, ventilation, or electrical systens 
to the various compartments. It, therefore, behooves the structural 
designer to make proper provision for these requirements without adversely 


affecting the structure. 


Naturally, the circle, "Slime of its simplicity, was the first 
type of opening chosen. Kirsch‘) analysed the effect of a circular 
opening in a uniform tension field; Tuzi'2) performed the analysis for 
the circular opening in pure bending; Howland and Stevdanen extended 
both the Kirsch and Tuzi analyses for finite width of plate and intrceduced 


the analysis for complex bending. wot’ and Neuber ‘>? 


6 
for the elliptical opening in pure bending; Inglis’ i a similar analysis 


performed analyses 


-or the elliptical hole in a uniform tension field; Wolf also extended 

his analysis to cover the ellipse in complex bending. Quite recently 
Greenspan‘”) introduced re analysis for the flat ovaloid in a uniform 
tension field. In a forthcoming maplicheaen”? Joseph and Brock develop 

a theoretical solution for the flat ovaloid in pure bending. This thesis 
is devoted to the latter problem. The theoretical solution in the AP-ENDIZ 
follows the method proposed in (8). The effect of finite depth is evolved 
by an analysis of experimental data rather than by a rigorous mathematical 


method. 





* Mumbers in brackets refer to the Bibliography. 
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It is the purpose of: this thesis to correlate data from photo~- 
elastic experiments witn a systematic series of modsis with varying retios 
of length to depth of opening and depth of opening to depth of web, sub-= 
jected to pure bending, in which the openings are Slut ovaloidss i.e, 
openings cemposed of rectangles with semicircles mounted at tae long- 
itudinal ends. This correlation includes a check on the mathematical 


solutions developed in (2), (3), (8), and in the APPENDIX. 








PROCEDUE 


The models were prepared using a nominal thickness of 3/8 
inches clear Catalin. All models were cut from the same sheat of material, 
and a standard tensile specimen was inciuded. The tengile specimen was 
used only to determine the fringe constant of the particular sheet to 
ensure satisfactory sensitivity as well as to provide means for obtain- 
ing checks through the equations of static equilibrium, if necessary. 


These models were machined using high-speed steel end mills for all cutting. 


In an effort to minimize the time~edge effect, polymerization, 
a specimen was tested immediately after machining, allowing only surfic- 
fent time for machining stresses tc disaipate. A white light polariscope 
was used for this purpose. The specimen was fitted with end brackets 
and suspended in an A ~- frame, 23 shown in Fic. I, which subjected the 
entire specimen to pure bending. Suspension *as accomplisned by vertical 
links fitted at each end with cireular pins of drill stock to assure as 
nearly point sup»vort and loading as was feasible. The A ~ frame was loca- 
ted in the center of the circular polarisccne using a mercury vapcr lamp 
with a Wratten #62 Filter. The resulting light was very nearly mono- 
chromatic with a wave length cf 5461 A°. The resulting fringes were 


alternately light and dark giving excellent contrast. 


The specimen was first suspended with brackets and no load. 
The analyzer was rotated to give a dark field. A picture was taken in 
this condition, hereinafter referred to as "no load™". The negative 


was then developed to assure proper focus before load was applied. 


== 

















Load was then applied, and a new picture taken. The load was previously 
determined which would give a maximum stress of approximately 100 pai., 

well within the elastic range yet sufficiently bigh to give abcut seven (7) 
fringes. The analyser remained in the same position ss in the no load 
condition. The negative was then Geveloped to assure satisfactory rocus 

and contrast before the specimen was removed for further machining. By 
following this method, it was possible tc complete a run, five (5) specimens 
- all with the game ratio of length to depth of opening, but with varying 
ratios of depth of opening to depth of web, in about four (4) hours which 


minimized the time~edge effect. 
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RESULTS 


Tne results of the experimental data are incorporated in Figures 
II, III, IV, and ¥. Figure III is a plot of stress concentration factor 
at edge of hole versus the ratio of length to depth of opening for contours 
of depth of opening to depth of web. Figure IJ is a plot of stress con- 
centration factor at the edge of the hole versus the ratio of depth of 
hole to depth of web for contours of constant ratios of length to depth 
of opening. On this plot the various experimental data is shown. Simil~ 
arly Figure IV and ¥V show respectively the same plots with stress concen- 
tration factor at the edge of web as the substituted variable. Table I 
is e compilation of experimental data from which the aforementioned curves 


have been plotted, 








TABLE 


Summary of experimental data obtained photoelastically 


Ratio of length Ratio of depth Stress Concentration Factor 
to depth of of opening to 
opening depth of beam 
¢ be 
1 2,06 af 
eae. ! "ee: 9 Ff Meee, |) Pe 
1.0 a) a. 2 aT 
weer MMI OD ne OR 
--- 1 
0,100 1.430 1,045 
Q,200 1,440 1,016 
1.5 0,300 2520 1.245 


2-5 0.300 42300 142954. 
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DISCUSSION OF RESULTS 


The curves shown in Figures II and III, and IV and V respectively 
represent the cross-fairing of the data incorporated in fable I. The results 
of Figures II and III are consistent with the well-verified facts that 
a small crack in the longitudinal direction of a bexum causes a stress con- 
centration of unity and a small crack in the transverse direction of the 
beam causes an infinite stress concentration. This character is shown 
again in the slope of the curves of Figures III. It is not surprising, 
either, that the circular opening causes a greater stress concentration 
than any of the other flat ovaloids (no ovaloids transverse to the axis 
of the beam were tested). This result is similar in content to a compar~ 


ison of (1) with (2) and of (6) with (4). 


The curves of Figures IV and V are not so easily explained. 
At the outset there appears to be no rational reason for the apparent 
maximum stress concentration factcr at the edge cof the beam for a ratio 
of length to depth of opening of 1.5. The almost perfect agreement of 
the experimental resulta with the theoretical solution for the value of 
the ratio of depth of opening to depth of beum of zero and Sor ail values 
of such ratio for the circulzr opening lend confidence in the remaining 
experimental work. On this basis the authors feel safe in predicting that 
the maximum augmentation of stress caused by a small opening of this type 
will occur when the ratio of length to depth of opening is in the vicinity 
of 1.5. This same confidence enables the authors to conclude that the 
ratio of depth of opening to depth of beam should lie between 1/4 and 1/3 
if no material increase in the stress concentration at the edge of the 


hole is desired. 


ios 








CONCLUSICNS 


1. The stress concentration factor at the edge of the hole varies 
inversely with the ratio of length to depth of opening and directly with 


the ratio of depth of opening to depth of bean. 


2. The stress concentration fector at the edge of the beam varies 
diractly with the ratio of danth of opening to depth of beam, ‘ut in some 


complex menner with the ratio of length to depth of hole. 


3. The critical value of the ratio cf length to depth of opening 
falls in the vicinity of 1.5. The critical value cf the ratio of depth 
of opening to depth of beam lies in the neighborhood of 1/4 to 1/3. The 
critical value is that value at which a small increase in the aforementioned 
ratio ‘gs accompar:ied by a relatively lurze increase in the stress concen 


tration factor at the edge of the hole. 


4. The thecretical solution developed in (8) and in the APPENDIX 


is satisfactory for purposes of engineering design. 


5. The photoelastic nethod is well-suited for inexpensive and rapid 


determination of stress concentration factors. 








eS 


RECOMMENDATIONS 


l. Further experimental investigations for values of the ratio 
of length to depth of opening in the vicinity of 1.5 seem highly desir- 


able to fix this critical value. 


2e Further experimental and theoretical investigations in this 
nature to determine the effect of reinforcement of such cpenings and the 


nature of complex bending; that is, bending with shear, seem desirable. 


3. In the absence of a olearly defined "no load® picture in a photo= 
elastic investigation, it is recommended that small but appreciable stress 


be introduced to clarify the basic pattern. 


~15~ 

















APPENDIX A 
DETAILS OF PROCEDURE 


Description of Models 
Four (4) basic models were prepared. Each was fitted with an 


interior hole 1/4 inchea in the transverse direction, but the longitud~ 

fnal extent were respectively, 1/4 inches, 3/8 inches, 1/2 inches, and 

5/8 inches. The distance between suspension pcints was maintained at 11 1/8 
inches; the distance from suspension points to the loading points was 
maintained at 5 inches. Each of these models was successively diminished 

in depth in order to obtain the increasing values of the ratio of depth 

of opening to depth of beam. Only in the case of the circular opening 

was it necessary to recut the interior hole. The particulars of the verious 


medels and the loads to which they were subjected are summarized in Table II. 
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TABLE II 


§ a fx i ticulars 


-18- 


Hole Dimensio Depth of Be b/c Load a/b 
1/4" x 1/48 22240 0.1115 85.35 1.0 
1/2 =x 1/2 2.094 0.239 $5.35 1.0 
ie x 1/2 1.372 0.564 46.00 1.0 
Iy2 x 1/2 1.030 00485 22.00 Le 
1f/2_x }/2 0,833 0,600 be00 Lew 
1/4 x 3/8 2.500 0.2100 Soo ie 
1/4 x 3/3 1.2500 0.200 20.00 1.5 
1/4 x 3/8 0.833 0.360 9.00 ine 
1/4 =x 3/8 0.625 0.400 £20 1s 

3.50 Q 1.5 i 

1/4 xl/2 22500 0.100 85-35 Pao 
me =x i/2 1.250 0.200 20.00 2.0 
1/4 x2 0.833 0.300 9.00 20 
1/4, x1l/2 0.625 0.400 4.00 220 
ib rele ‘ 500 
V4 x 5/8 22500 0.100 56.90 2.5 
1/4 x 5/8 1.250 0.200 18.50 205 
1/4 «x 5/8 0.333 0.300 9.50 Bins 
1/i x 5/8 0.625 0.400 6.00 os 

8 0.509 9.500 5200 225 





Kethod f alys 


Standard photoelastic methods were used, supplemented where 
required with direct analytic methods. A transverse section was chosen 
a minimum of two lengths removed from the longitudinal end of the hole. 
The condition of stress at this section should, according to St. Venant 
be removed from the influence of the hole. This secticn was designated, 
accordingly, "pure bending", see Figs. VI, WII, VIII, and IX, and was the 


reference condition for determining the stress concentration factors. 


Examination of the photographs, Figs. X, XI, XII, and XIII, 
revealed that the maximum stress at the edge of the hole occured at its 
vertical uxis of symmetry, designated "diameter"; whereas the maximum 
stress at Sua edge of the beam occured vertically above the longitudinal 
end of the opening, designated "end of hole.* It was therefore necessary 
to determine tne fringe order existing at tne free boundaries in these 
two locaticns. Since a frings is the locus of points of constant differ- 
ence between principal stresses, one of which must be sero at a free boundary, 
no other experimental information was necessary for the evaluation of actual 
stress existing at these free boundaries. Fortunately, the "exact," math- 
ematicel, stress analysis for pure bending also raveals that one principal 
stress {3 also zero everywhere, which fact obviated any additional inform- 
ation for the pure bending ccndition as well. Thus the desired stress 
concsntraticn factor was merely the ratio of the fringe order actually 
mensured at a given point to the fringe order determined at the corres~ 


ponding depth in the pure bending condition. Both frings orders referred 
t 





+ Ref, (9), pe 31 
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to above must be the actual valuc cetermined by the photograph corrected 


for any residual stress in the "no load® condition. 


This correction is show in Figures VI, VII, VIII, and IX as 
the difference between the dotted and solid lines. The case of pure 
bending was obtained by plotting integral fringe orders, solid lines in 
the photographs, at the proper fraction of depth from the bottom. The 
upper and lower halves of the beam were handled separately, since, in pure 
bending, one half would be in tension and the other in compression, with 
zero stress at the middie surface. The best straight line, guided here 
by the theoretical linear diatritmtion of bending stresses, was passed 
through the experimental points in each half of the beam. It will be 
observed in Figures VI, VII, VIII, and IX that in no case did these lines 
pass through zero at the middle surface. Shifts of both halves of these 
lines affected to make the lines pass through thia zero point and to inter 
cept the same value at the edges of the beam. It wiil be noted from the 
figures that in most cases this was almost a parallel shift of lines with 
approximately the same values for both halves of the beam. This amount 
of shift, in all but the shallowest mcdels, was of the same character as 
the load stress and less than fringe order 1. This represents the residual 
stress in the specimen, which from the density of blackness in the "no 
load® picture, must be very nearly constant over the section. It may 
be noted, however, that in the shallowest —.: is appreciable resi- 
dual stress, more than a single fringe order and of the same character 
to the load stress. This condition is caused by the bending effect of 


the end brackets and admits a more straightforward correction procedure. 








The Giamoter® and "end of hole" curves have this same correc- 
tion factor applied to them in all cases except for the circular ovening, 
Figure VI and X. This procedure was followed since the dengity of 
blackness on the “no load® picture appeared to be the seme throughout 
these specimens except for insignificant surface blemishes. The circular 
openings retained machining stresses for a considerable period of time. 

The fringes representing these sachining stresses were concentric rings 
which permitted a plot of fringe order versus frection of depth from the 
bottom, see Figure VI. These values were then added algebraically to the 
experimentally determined values, previously plotted. The character 

of the machining stress was easily determined oy noting the *"bulgine® 

or *depressing® effect on the fringes in the loaded condition. For example 
when fringes observed in the compression half of the beam were "depressed® 
toward the hole and bent away from the diameter toward the end of the hole, 
thus increasing the fringe order at the diameter of the hole, the machin- 
ing stress was obviously compression and should be subdtrected from the 
experimentally determined value to obtain the actual fringe order caused 
by the bending alone. Conversely, if, in the compression half of the 

beam, fringes were *bulged® away from the edge of the hole, causing dips 
on either side of the dianeter, the machining stress was tension and 

should be added to the decor techibhl ie determined value to obtain the actual 


fringe order caused by the bending alone. 


A polar plot to obtain the experimental fringe order at the 
diameter of the circular hole wae used in an effort to ascertain this 


peak value with greater accuracy. This is shown in Figure VI. The 








photograph was tacked to a drawing board and the center of the hole deter- 
mined. Radial lines were then passed from the center of the circle through 
the intercepts of the fringes at the circumferance. The angles were then 
transferred to the plet and the correct value of fringe order was plotted 
radially on this iine. Finally, a smooth curve was passed through the 
points and the value of the fringe order dete-mined on the diameter ex- 


tended. This gave the desired end point of the curve. 


~Lan- 
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APPENDIX B 


Sample Calculations 


The data are Summarized in antirety in Table I of RESULTS. 
Figures VI, VII, VIII, and IX, referred to in APPENDIX A., are shown to 
4{llustrate the methods of analysis and typical calculations of stress con- 
centration factors. Ons sample from each series of models is included; 
that is, one for each retio of length to depth of opening and each ratio 


of depth of opening to depth of beam. 
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APPENDIX C 


Original Data 


The original data is presented in the form of photographs of 
the *®no load® and load conditions. Each page of photographs includes all 
the pictures taken for a given series of modelss Figure X for the ratio 
of length to depth of opening of 1, Figure XI for the ratio of 1.5, Figure 


XII for 2, and Figure XIII for 2.5. 











b/C = 0.600 No Load b/c = 0.600 Load 





b/C = 0.364 No Load 





b/C = 0.36" Load 


ne 


b/G = 0.239 No Load 0D Sz 
b/G = 0.239 Load 





t 
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rr SS 
——— 
ae 
s/o = 1 ; - aati 
\ 
\ a 
FIG. xX 7. 
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b/C = 0.1115 Los 





-~3i- 








Pe 
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b/G = 0,4 No Load b/C = 0.4 Load 





b/C = 0.3 No Load p/C = 0.3 Load 





b/G = 0.2 No Load b/c = 0.2 Losd 


4 on 


ots i, 





b/C = 0.1 No Load a/d = 1.5 b/C = 0.1 Load 
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b/C = 0.5 No Load b/S = 0.5 Load 





b/C = 0.4 No Load b/c = 0.4 Load 






b/C = 0.3 No Load 


b/c = 0.2 No Load 


a 











b/F = 0.5 No Load 





vb/C = 0.5 Load 





b/C = 0.4 No Load 


a/b = 225 


FIG. XM 
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v/C = 0.4 Load 


PRE: oe 





v/6 = 0.3 Load 





b/G = 0.2 Load 





b/G = 0.1 Load 

















FIG. XITv 


X= CONSTANT KYA \K 
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Nomenclature 


(x,y) 
Neo eet 


z= 


w= 


G 


h= 


Ae 
HEH 
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<6 
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ee 
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A 
A 
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Ans Ens @n» Dy 


Hi 


Jo 


APPENDIX D 


Cartesian co-ordinates 

Orthogonal curvilineer coordinates 

x Ly 

a + ip 

ev 

stretch ratio 

et P 

angle in z-plane between tangent to curve 
Bp = constant and x-axis 

real parameters 

oo ob - 6 

potential functions of complex variable z 
potential function of complex variable 
Airy stress function 

normal strefeses; 0,18 normal to surface for which 
x =constant, etc. 

shear stresses 

comolex constants 

apolied bending moment 

half-height of opening at x =0 

B- tf? Rk 

moment of inertia of beam 

bar above, indicates conjugate complex 


ctress concentration factor 


hae te 








Levelopment of Theoreticzd Solution 
Consider a function: a= fete) (1) 
which transforms continuously points in the Z{a, 8) pune 


into points in the 2(x,y) plane. Here z= x ¢ ly and 


rn oe Me 


ae 


= 


In purticular, see Figes. X¥(a) and (b), points on tne 
circle of radius ae with a hele constunt, cefinead by a 
mosition angle gp , will transform into voints (x,y) ona 
closed curve defined by the functional form (1). The reéelal 
line defined by Bp = constsant trunsforms into a curve orthogonal 
to those for which a= constant. The case of @ = O if in- 
portent; it defines the unit cirele, y , in the @& -nlane 
and transforms, by the proper selection of the form of (1), 
into the shape of the openine to ve studied. A special symbol 

» will be assigned to those C-valuesc which lie on this 
Circle: 


A = ptf (2) 


M@erthnat the conjugate of A:A= ew Bos 1/7 a 


The transformation defined by Ecuation (1) describes the 
relation betiveen an increment dz in the z-plane and an increment 


dv in the y-plane. The ratio, exprersea in polur form: 


2 ne ae (3) 


id 
Ty 


2 =. thie cdevelovoment followe that of (8) and wes lurcely guided 
Peer. dg. 5S. Prock. 











defines hk the streteh ritio, and X , the ancle betreen the 


tangent to the curve fio Ce a and the x-axis; see [lol , 


Article 2.32. han’ Y are obtained by differentiating 


adi q a d, & 


Pauatins (3) and (4); 


Gs ay (h. az 









where the prime indlcates @ifferentiation with reepect to 1 
variable in the purenthesis. 


yy 


Thus, the conjugate of (5) is: 


h en ix 


povend 


7a ae) (67 


G c 
provided (1) hee real coefficients and parameters. 


Wultiniication of (5) and (6) gives: 


He 28 ( Zee) Fan 2 eee ee (7) 


oe 


The division of (5) by (6) fives: 


9 clk = zife) yee P (8) 
Ze) 


In general, in plene problems of elusticity, the teen 
miaue which follows can be anolied to any renl rational m 
Meme tunction 2 =f(¢); thut is, where 2 18 the r.tio ope 


polynonl..te in G. maivin peal coerticiente. 


- 39 = 








Consicer the mapcine function due to Ureens oun 7 


ry 


4 3, (9) 


mnere Ss, t, .nd fi tre rewl constants thet are Geter. 


a 


BIG) Ae 


Mfc 


from the proportions of the onvening, and z = x 4¢ ty and 


a = a+i§ 


Then: «x 4 ty = set. ef B 4 tema . e71B YZ RemFa . em- 4B 
= (se@ £ te~2)caor px Re “Gcoc3. f i(seS- te~2)sin 3 
~ 1Re7 a sin3 p. 


Ecuating the reai and imaginary narts separately: 


“FP 
ty 


(sa@ ~ te) CO6B Z Re7 2a cos 8 
(10) 


Yee (ce ~ te?) etn, = Re77@ gn} B 


% 


when g = O, Pouations (10) define the general ovaloid opening: 


x= (s # t) cosg # R ens} 2 
(19a) 
y=(s -t) sing - R ein? g 


beveral special Gases are of interest: 


(a) Wnen R =O, an ellipse is cefined. 


ee 
HW 


(cb) when 4 =O, a circle of radius @ 1s define 


we) Wien O and F.< 0, a scuare with Pillets at the 


c 
Hi 


corners is avvroximated. 
{a} When t = 0 and RF C, the Snouare is rotated cee 


in the general case, by the proper choice of s, t, and R, the 








ecuations can Le meade to anproxima te {ne Co. vonec trict 
openings, evecificaily tne one Gomeesea Gt 2 rectsne.= = ita 
eemicircler at exer en. 

2 iin 


From (9): 2) (Ge oe ee ee) 


: a 
oo 
Now hand YW can oe determinec from (7) and (8). When 
a= 0, corresponding to tne opening, the suveerint zero 


me used enc from (7): 


=49( A)z'( 4) = (eg = 2p a - - <7) 


: ws Ni 


ofr 


2 z | 
=e + t° FOR” ~ (26% - 6Ht)cos2 8B - SAscordP 


In order t> obtain strernses in nlaine problems of 


elusticity by uSings the complex variable method, tvo votential 


fFunetions muss be deterained. These funetions can be eaxoressed 
eo: 
x n i) n 
Sz) = > Aaa ee 
{ n § cas 


20 


x ao 
Y(2) ra Ron of Ze Sey 


_ 


where, Ans ns Py, and Db, are constunts, posFfitly compilers 
“ne solution is thus recducec to tne evaluation of there con- 
eypants. the stress confition at one bouniry, infinity eee 


cetermiine A, and 5 the strercc concition at the other bouncary, 


"2. 
——— 


the openini, will cetermine al and bl. Onse the tyro potential 


aaah | 1 

functions sre completely ceterminecd, the stresses tn the onlane 
Je v He 

may be determined py diffrerenti..tion. Thus, there functions 


2re countern.rts of the Alry etre:s function. In fact, the 


or yee 








Siry strers function corres oniing to tnese potenti.l functions 





is given by: 


Pix, yy) = ee ee: of | 3 (2) r fVerea] (144) 

The stressere cun ve found in ters of the two functions: 
= 2? , WA ‘ 2)| 
o,, 4 o.. [ 4 (2) £ Pz) 

(15) 


ay - ox  2tyy 212 Gra) x Pz) 


| z Hy 
The usual transformation equatione ive: 








“Ss, - o. + cit je = Oj Oo mae elt, (16) 


By combining Scuatione (3), (15), and (16): 
o, +0, * 2LO')¢ 6+ (2)| (17) 


-o, + 2ha,,= 2 [2(2) Or(2) + Pi(2)} 2442 e218 


“8 (2) 


where the functions of 2 are to be exvnressed in terme of 


mec ,p ): 
Piz) = Pkla)l= 9 lo 
Pos) = Pei z)| = o (x) 


Hifferentiz:tion is with reapeat to the varluble in Surentae- 


(18) 


secs 


hs Ref. (10) , ». 135 


- lo = 








WG eae (19) 


>" (7) = gp" (S} : _ Q's e G 
Bora [2a 


if, at the opening, the normal stress on anc the she.r 


5) 
+ 


a9) 


srese Liat are zero, which is the c se at a free boundary, 


1% can be shorn that: 


gp (A) + ~-Sa54- o (2 )6 eae (255 
r N 


Tie functions 9 (2) and W ( Z) may be evaluated 


with the use of tne Cauchy Theorens regurding poles and resi-~ 


dues.? 


Prom (12), these are: 


Theoren T. If f( Z%) is continuous in the closed region 
|x| = 1 ané analytic in the interior «1th the possible exce = 
tion of the »voint Z = 0, where f( Z%) hus the form 

“~ ee fi %, 

ge) = ee t AA circa) 


ana «{ 2%) is analytic: then: 


where Y is the unit circle y = el B 





Theorem II. Li f£(Z ) if continuous in the closed region 


\y\2 1 and analytic in the rerion exterior -o A” , with 


the nossibdle excention of the vnoint y= 


has the forn 


oO , where fOr) 


| By 
f ( S } = Ay po AGG vs As c a oo @ a Le me 7k 
ihens 
1 f(A) ie 
eo [ia = =o) A. tA, St a fe scam ceeds 
POR: | ale 
Consider Fauations (13): 
a = ao na 
@'(2) 5 a DAZ ae a _ “at 
, - n-2 = n(n+lja, 
« 
ze aroll = nb, 
i a a 


Now let 2 aporoach infinity in the z-direction: 


o>! (z ae fe | > nie | cn 
§: (2) #2 | z hao 
Oi(z) a [= naye-2| aan 


The boundary conditions at X = Pfare: 


Pe oe me 


L= eo 


(219) 





Inserting E-uations (21) into the first of =qaustions (15) 


along witn these bouncdury conditions gives; 
= n=1 > n 
= ie 2| 2 na n(x x iy) Z Zona (x - iy) bees 


This reculres th.t Ay, =0 fer all values of n excent ¢. For 


2 
2 ey =2 |2iagy - 2ihay| 
_a kf 
o2- Siias 


Inserting Pceuations (21) into the second of Eruations (15) 
glves: 


© 
4 e 
) 


“y= 2 lx ~ Ly)+2+1-An¢ a nb, (x # al 


ae = &O 


“his recuires that Bn =Q far wll values of n extent 27 ee 


Ey = une ZABRY 
Fy 2 ire - Z1Boy | 


7 a. 2 
BF a 1 ayy the 


Therefore, for nure bending, regardlesa of the shape of tne 


opening, the potential functions are: 
b(z) 


y es > 
Pz) Bt 


On the boundary of tnie opening: 


1) 
ti 
= 
tO 
“WN 
a lpais 
1) 
es 
™ 
Re 
pm 


i! 
i 
t 
eet 
t 
~N 
M68 
of 
ee} 
Re 
Bs 
S 


aoa 








mp (A) = te ee +8a)* + Z. eee 
eu! A 2X | (8 444 +H)" 
A 2X 
rd ier. = ni I 
~~ (= +5 \ U4 eee 
4h t RR \é ee by 
y (A) =- ==(s At = +~x ) - 2. Pps & See oat 
5, * 
1 mb L a Z = a 
% ae Bea ee AS ee a aw ginal er ta slr = 


The boundary condition (20) besomes: 


4 Zz as 
A+ tA” + RP Ads = na 
yay EMER A(R ve ne Pn pe tlnenegae 
a 40 A (= +tA + RAW) 
: ran) oO 5 
oe eae +tr + FA Bae = rae ge, oo oe oe ow ee ow oe = O 
8I r 8 3\n 
(= +t + ae 
(egy 
Expanding and collecting terms; 
t . £ 2) e 
vw (A) =- “7 (8. ee Ais -st) - a= (28t 42Rt = t” - ZRs) xX 
i 4 4162.6 ‘SH 26% +28t See 
~ =-(2He - 2Rt) A +=2-R7° QA - = (e------+3----- ) 
St eee SI SI 2 
9S 8 Ee 
31 A (Seta hx) 3 
7 sa oom coms So ee Dy on om an an eae =) 
+ g ome 6! 
Te a) 
r 
Me) 2 =0, of x) hus the for: 
‘se iu 2 aes Zs a 
is rt -s ie c2RS | Ra 
Cy) a ety) gh(2ee tert ss? 4a, 28s, 2 
ye) = Eg) +o - Ol xy” Be 








rrom Theorem f[: 


i @(A )aA ih, 2att 2Rt -s* if 2Ps Fe, 
wt of T8225 ~ sa(sShech nh) - ah. ook. G1 
2 mh hax BT = BI x2 ee 
A 
. 
for ly | > 1, 

Bs _f 39(A)0).. _ 1h cats, See 
On 1 2 BI x 31 a 

y for |Z1I> 1 


At ~2~= eo , go (x ) has the form: 


ot Sol ae 
i e , \ 
gf Z) = a= 8% wo +d + 
BT ¢ a 
srom Tneorem IT: 
eo. Oech = 6 ee for |zl>l 
LT 4 AW-Z% Sy 
y ADLA JAA —" 2 
a f “ (AEA = - %@ (Z% ) + in Bre oe for \ ae 
2 ai AWKZ OT 
Y 


From the first of the expressions derived fron euch theoren: 


iT .2 i ORy =e° , ce RE 
o( Zz) = a. ‘eG re ae + == + ae (2) 
5 S eee, 
rrom the second of the exnressione derived from each theorem: 
ii4[ 2 (32 2et+ 2Rt -s*, 2hs Pe 
vo (zg) =the So See + <- (25) 
BI y ae 


Comparing (24) and (25): 








Putting T =t # a - z0 
= its g% a ay a 
9 ( a) Tt [ 3° a ae (26) 


But the conjugete of Ecuation (2°) 1a: 


= (dt) + ZAR gr(r) + WEA) = 0 
ys z' (A) 


Using taquation (26), this becomes: 


g ae 
. 18 | =f +7 AS +p 0 + 23 ee i +t oe : 
4t k2x BY 3 AD tg = oS ee 
A ee 
+ y fA) = 0 
. ec) 27 = HR 
“ae . +PA° ¢ oe as ta a D2-(5 4% 7° 4k Ne) 
4E "2k SA = ta” <3R 
At y = 0, y( &) nas the forn: 
- LE(_8 
y iy )= nT ‘372! + fe) 
From Theorem I; 
WoA)a is, s : 
eee tor Jel> 2 
25 1 “% i 2 


y 
oe. = - ,¥ (x ) has the forn: 
a 
iNe Be ee 


LR 
AT SZ ee - 38 


From Theorem il: 


4R 
; 3% - 2T - « 
da a & 
a WAIdAL -W¥ ( Z)+ eet tr + i y - ee 
eal AWZ 4T sx = tz een 
Y (oat ae’) | for | ame 


_ er. 








Therefore: 


HEE WR 
is] _s : b sf -2. - <6 — 
ter g atc see 


Using the first of Ecuations (17) und notin that at a free 


fof) 8G). 
z' (x) aes 


boundery a = OC ae 


(op ) =2($'G) b(2)| = 


) 
a=0O a=0 


N) 


But 


? Zz 11 = a ky 
af PN re a fom 7 a ae 
gl A) = zs . “2+ | and @ (¢)= - ae | =£2+ ie ar'| 


: 6 “~? t : c 
ot (Aj= AUS [= - £4 - “ji and @ '(4)= - ius| ~2tx3-un} 
47 Ae N uS NT 
ae Z 
g(A) SE = Se = or and 2'(t) = s-- t x%~ 3R n° 
Thus $ 
: 2T &R E37) oe 5 
to. ) <2 | ee. Se. 
B a= + a 2 
8 = “ - “kh 41 Eo te 3RX" 
x » 
(o,) = - 22 | Set SRec gts. 228! Jotne- (ent stse® detnspemennge| 
Pa=0 T Us +e tOR ) +(6Rt - Z5t)cens2B - Ghscos4s P 
(Zon 
(04) is a maximum when p = 90°. 
a=0 
Therefore: 
(o,) fi 3 7; 
ae 25° —-3Rs-1Z2R°+2t(4Re - 
CCR = 2-929. s 2 (2R:) 








If ‘sa. nou suvstitute Fauavtone (i) (97 G  eee 
(11) into Evvation (14), interrate me camiine mens ae 


get the Airy streee function: 


WY = ;° {oe act; ry wie ~ a =a 
rla,B) — ae { [ (02. PaOtt one lee +(32f - = -2 Tt -4R“)e 


Z 
e¢r Of oO OG at cam Tae: p= -¢ 
, (cae gee 2a seal “0 eee 
3 6 3 3 


to 


: Z a -5a 
st a -5a "ud pe "2 Bg 
-— --e + hte ‘| singe +[ R26 - --@ = pg ] Bing at 


her 








- 
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*valuation of Parameters 

It is necessary to ceternine the vniues of RAR, s, und t 
that will eneure that the mapping function (7) will suit thre 
actual opening when de C. Several conditions may be ing 
posed on (194) in orcer to cetermine Ff, 6, an. t such as: 
{a} sx tA Re 


(pb) ¢ - t # Re = vertical dimension = +b 


afte 
“% 


horizontal dimeneion = 3a 


$ 


(c) Srea of opening = (a - t&)b # pq = Jyax 

(dQ) Radius of curvature at end = ib 

Fauatlons xusscolatine the desired varlabler with the 
geometric properties of the openin. must no be ceveloprd to 
provide a systematic vrocedure for ovalue:tion. To simplify 
this oroblem we introduce the Greenspun notation: 


#(P 4) = a 


4 


(30) 
(P -}- &) 2 Ft 
or: 
sx t= P 
(39a) 
Ss = c = i 
Then, Equations (10a) become: 
X= P cose - R coe 3g 
(1. Ov) 
Imposing the concition of edcusl areas, ve obtain: 
eu 
fy vi y axis =) pes kp tor 2° 1S ( 208) 


Sapstituting (1Cb}) into (231), integrating, and @@imolityes 


FiveEc;: 





=(a - b)b # tb“ (32) 


a 


Reference to any text on differential calculus showe that 
the radius of curv..ture ut any point is given by, toe eee 


oregrion: 


we 
Radius = oe ----- (233 
y 
7 @ 7 4 —= ay — cy e Cp » it — dy! — Gy! * dB 
ee YY Fay da’ dx SM Y= ae ae To ag 
But, Ce =Q cosp - 3R cos 38 
and a3 = ee oe wk (ar aoe ek ee 
Ce P sing # 3R ain 38 
Hence: y! = -— & COS fp -_ BR GOs 3p ( 344) 


P sinks 2% ein 36 


Ba 2 ~ (* 5 
and y" = 4 Pot 27 -3BR(P Yc dooe 4 BY 6R(P-o)cos_ 2 P 


(P sing? 3R sin 38 )e 


(327 


Substituting (24) and (35) into (23) and simplifying gives: 


” 


Mle Pe jae, re P [ea 
ndius - 2bPoeol elo PoP BP Re Bone He" eos? PBR P+ Joost BOR” 


nov p= O for the lonzitucinui end of tne opening. 


Pmerefore: 


Padius = {4 = 3R0° = ib (37) 
p 4 9R 
Tne remaining two equatione are simply: 
Pos Ra (38) 
(, f R = bd (39) 








Considering only Mewuetions (32), (47 502) ee ee 
ae Bys8ten of enuations for the determination of PF, G, ami &, 
we nox chooge b= ? and involve the ratio of lenrth ta cepoth 
of openin:., u/b . The cystem of ecuatione become: 

Pf#£R =za/b 

eo 


4 be (740) 
ee fe = (a/b -1)/1 
i ee 
fc en =e 
P # GR 


Any system of more equations than there sre unknotns is ine 
Geterminate. ‘“'herefore, the four possible comuinstions of 
(40) taken three at a time vere solved for P, 4%, ant R for 
the particular case of a/b = 2. This particular case was 
chosen since it had veen used by Creenspan [7] as an il- 
lustration and a ready comparison wes available. The outline 
of the opening as defined by there four solutions are plotted 
in fig. AVI ang comoarea with the sctual opening for which 
the fit was cesired. 

Because the ratio of denth of opening to cevtn of beam 
is an important criterion in the determination of tne stress 
concentration factor, the deoth of openings seemed the mort 
logical first choice to be -eld constant. Next in invortance 
as governing variable was the length of opening. To complete 
the system of eoustions for solution for P, @, ane F either 
tne radius of curvature at the longitudinal end or tie ares 


must be chosen. A comparison of the fit ehnoved that either 


~ 53 - 








Si deevasnasvedvadsaenedvaussved seers uassserat fvetaesaaz seeessveesaedastaedseeaetieeveeraetisza 
EGGERS SERA ORES PEON SO OOOSONoo Boao Senate eater eee Ea 


sdggaawsas sssesseueecesscccc555555 erasers tere eee 


BESS CEES BASS Pee SSeS SLSR OBeaow BS CSEES CESSES ESE SSSR BESS GSS G Ase ] 
ele Suzi case eeasseaud uses aaseetetasteeaa factinaiii 
sam Seen Tit th ae BESTS Pees Seaes Se seeeense 206 stesee 
ear iittt SSS sestasesaseaazeet eee tgee teat szeed ete tot tezvazee 
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wroule be satiefactory, ond u comparicon of the two values of 
strese concentration fuctor emanined wit): Greensvan's value 
&s the vurdstick anoved thit either ty mid be satisfactory. 
A go0d cise could be made Tor noleing Gee sence erec eee 
Curvecure ecusl 0) tne lhuidiyv weistt of opening since thie 
system cuve values of stress concentration factor thet were 
Slightly high - pessimistic. When plotted on a curve eimilar 
wo Fic. Lil, however, the closeness of fit of values culou- 
latec by holding the area constant vas stizchtly superior to 
those valuee celculated by holdinc the end radius of curvature 
ecual to the half height of openine. The s.se or calcr Tava 
the valuee of PP, %, and R when holdings the orea constant wee 
a further consideration. 

Consecuently a cysetenm of equations consisting of tne 


firet three of Feuations (4::} vas chosen: 


P= a/b - R 


Gs 1 = (40) 
Pij= IRE = 2 (aft -1i)/41 


From which: 


R= - -(s/b41) ¥ 2 (46a ) 





(a/bdl <= Bie) aye 


hese valuee of P, G, und RB are vlotted versus a/b in Fig- 
AVIil. 

this should not be construed to mean that these values are 
precise; they cre merely the best choice thet the autnore could 
make with existing information. tome cther geometric property, 


euch «us the perinheral Lleneth, may warrant consider:tion. 
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